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ABSTRACT 



Oh 



Context. Spatially resolved brown dwarf binaries provide the unique opportunity to determine the orbital motion and dynamical mass of the 
system independently of theoretical models. This is a key ingredient for the calibration of the mass-luminosity relations and evolutionary 
Q tracks. 

^ Aims. We have monitored the benchmark L dwarf binary Kelu-1 AB over the past 3 years to derive improved spectral types and luminosities 
C/^ for the individual components. The astrometric measurements enable us to compute the orbital parameters and result in the first dynamical 
, ^ , mass estimate for the system. 

Methods. We obtained resolved high angular resolution, near-IR images with HST and the VLT/NACO adaptive optics instrument in the J, 
H and Ks bands. In addition we used NACO to achieve the first spatially resolved mid-resolution spectra in the H- and K-band for a precise 
^ spectral type determination. 

\^ Results. The near-IR spectrum of Kelu-1 A reveals a distinct dip in the H-band providing evidence that Kelu-1 A itself is a spectroscopic 
binary. We derive spectral types of L0.5 + 0.5 for Kelu-1 Aa and T7.5 + 1 for Kelu-1 Ab. Kelu-1 B is classified as spectral type L3pec +1.5. 
From the relative orbit, we derive an orbital period of 38^^ years and a semi-major axis of 6.4^^3 AU. This yields the first dynamically 
determined total system mass of 1 77^53^ Mjup for the Kelu-1 AB system, with the uncertainties mainly attributed to the error of the system 

^-H distance measurement and the yet missing apastron passage. The derived properties of Kelu-1 AB allow us to test common theoretical models. 
The comparison of our results with color-magnitude diagrams based on evolutionary models yields a slightly revised age estimate (0.3 - 

00 0.5 Gyr) and a discrepancy between dynamically and theoretically derived masses, stressing the importance for further dynamical mass 
determinations of brown dwarf binaries. 
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1. Introduction 

Kelu-1 was discovered as one of the first single free-floating 



brown dwarfs (BDs) in the solar neighborhood by Ruiz et aL] 
(|1997 ). Originally classified with a spectral type L3 in the near- 
IR (Knapp et al. 2004 ) at a distance of 18.7 ± 0.7 pc and with a 
large proper motion of 285.0 ±1.0 mas/yr ( |Dahn et al.|2002 ), 
it has been intensively studied since. More than 100 refereed 
publications listed in the SIMBAD database made it a kind of 
a benchmark object for L dwarfs, especially in spectral classi- 
fication schemes. Particular interest was devoted to its appar- 



ent overluminosity compared to objects of similar spectral type 
([Martin et al.||1999b[ [Leggett et aL||2001[ [Golimowski eFaT 



2004a), and to its photometric and spectroscopic variability 
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with a 1.8 h period in ( [Clarke et al.|2002|[2003] ). Since the 
presence of the Lii absorption feature at 6708 A confirmed 
its brown dwarf nature ( [Ruiz et aL 1997t [Kirkpatrick et al. 
[1999| ), the brightness could only be explained by a surprisingly 
young age (-10 Myr, Golimowski et a l. 2004b ) or unresolved 
binarity, as first proposed by Martin et al.^ ( ,1999bj ). Due to the 
free-floating nature of Kelu-1, its age is difi&cult to assess. In 
a first estimate [Basri et aL ( 1998) constrained the age on the 
basis of its Lii absorption strength to 0.3-1 Gyr, which was 
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later slightly revised by |Liu & Leggett] ( |2QQ5"] ) to 0.3-0.8 Gyr. 
This meant that a very young age was unlikely the cause of 
the overluminosity. However, Hubble Space Telescope (HST) 
NICMOS observations of [Martin, Brandner, & Basri| ( [l999a| ) 
could not resolve any companion in a first attempt. 

Six years later the mystery was solved. A binary compan- 



ion was detected independently in March ( Gelino etal.^200 6 ) as 
well as in May 2005 ( [Liu & Leggett|2005 ) with the Keck Laser 
Guide Star Adaptive Optics (LGS AO) system. New (photo- 
metric) spectral types of L1.5-L3 for Kelu-1 A and L3-L4.5 for 
Kelu-1 B were determined, respectively. We were able to con- 
firm the binary nature with our own HST/NICMOS data ob- 
tained in July 2005 and could already detect an orbital motion 
with an increase in separation by 0.015'' and in position angle 
(PA) by 0.43° in five months. This made it a new promising 
L -dwarf binary target for dynamical mass determination, due 



to the first predicted reasonably short period of ~ 35 years ( |Liu 
[&Leggett'20'0T). 

Observational studies of brown dwarfs play a crucial role in 
the verification of evolutionary models and the theories of in- 
ner structure and atmospheres of substellar objects. As brown 
dwarfs never stabilize themselves on the hydrogen main se- 
quence, there remains an ambiguity between the temperature 
or luminosity of a brown dwarf and its age or mass throughout 
its lifetime. However, spatially resolved binaries are expected 
to be coeval, thus removing part of the degeneracy in the mass- 
luminosity (age-temperature) relation which in turn can help to 
interpret their physical properties and test evolutionary models. 
In addition they provide the unique opportunity to determine 
the orbital motion and therefrom derive the dynamical mass of 
the binary and its components independently of any theoreti- 
cal model. Although these dynamical mass measurements are 
strongly needed for very low mass stars (VLMS) and brown 
dwarfs, only few such measurements have been achieved up to 
now and large uncertaint ies remain (see|Zapatero Osorio et al 



2004; Bouy et al. 2004; Brandner et al. 2004; Stassun e£al 
2006; Ireland et al. 2008; Liu et al. 2008; Dupuy et al. 2008 



Therefore we started a detailed ground-based high- 
resolution adaptive optics photometric and spectroscopic ob- 
serving campaign of the Kelu-1 AB system, including a long- 
term monitoring program to measure the orbital motion and 
derive a meaningful determination of the orbital parameters. In 
this paper we present the first dynamical mass estimate for the 
binary system Kelu-1 AB independent of theoretical models. In 
addition we show the first resolved spectra for each component, 
leading to improved spectral types and evidence for a possible 
unresolved third component in the system. 



2. Observations and data reduction 

The results of this paper are based on a variety of observations 
with diff'erent telescopes. For infrared photometry and astrome- 
try we used our own HST/NICMOS and Very Large Telescope 
(VLT) NACO observations as well as archived SPITZER ob- 
servations. The infrared spectroscopy was carried out by means 
of own NACO data. Each observation will be consecutively de- 
scribed. 



2.1. HST/NICMOS Photometry 

KELU-1 AB was observed on July 31, 2005 as part of our 
HST spectral diff'erential imaging program GO 10208 (PI: W. 
Brandner), targeting 12 isolated L dwarfs with no known com- 
panions so far. The sources were observed with the NICMOS 1 
(NICl) camera, providing a high-resolution pixel scale of 
0r0432 and a field of view (FoV) of ir'x ir\ in the two 
narrowband filters F108N (1.07 yum) and F113N (1.12 yum). We 
received 4 diff'erent images for the object in each filter. Two im- 
ages were obtained in the same orbit, but at diff'erent positions 
on the detector, to be able to optimize PSF sampling, reject bad 
pixel and gain redundancy against cosmic ray events. All data 
were acquired in MULTIACCUM mode and the total integra- 
tion times per filters were 2560 s (F108N) and 2816 s (Fl 13N), 
respectively. 

The HST data analysis of the resolved binary is based on 
pipeline reduced frames as provided by the HST archive. To 
derive the magnitudes from the HST/NICl images we used 
aperture photometry with the IRAF phot routine in the daophot 
package. With a separation of 0.3'' the distance between our 
components was less than the defined 0.5'' ( = 1 1.5 pixels) sep- 
aration for well-isolated sources. Thus we used an aperture 
size of 3 pixels and corrected to 11.5 pixel photometry using 
an aperture scaling factor derived from TinyTin|^ PSF mea- 
surements. Then we adjusted our results to a nominal infinite 
aperture and converted the count rates into flux using the most 
recent photometric keyword-value as provided by the STScI 
webpag^ For better comparison, we finally converted the flux 
into the Vega magnitude scale, using the flux zero points of 
1937.0 and 1820.9 Jy for the F108N and Fl 13N filters, respec- 
tively. 

We also reexamined the first HST/NICMOS observation of 
Kelu-1 from 14 August 1998 (GO 7952, PI E.Martm). We re- 
trieved 2 images from the HST archive which were acquired 
in MULTIACCUM mode with a total integration time of 448s 
each. 



2.2. VLT/NACO 
2.2.1. Photometry 

Follow-up imaging observations were obtained with the adap- 
tive optics system NACO at the ESO VLT/UT4 (Yepun) in the 
/ (1.25 yum), H (1.65 yum) and (2.15 yum) broad-band fil- 
ters on April 28, 2006 as wefl as in L (3.8 yum) and NB4.05 
(4.051 Jim) band filters on May 26, 2006. Between May 2006 
and June 2008 we acquired 6 additional epochs of Ks broad- 
band images within the framework of our ongoing astrometric 
monitoring program, as listed in Table[T] 

We used the CONICA S27 (for JMKs) and L27 (L, 
NB4.05) camera, providing a FoV of 28'' x 28'' with a pixel 
scale of 0'.'0271. Wavefront sensing was performed on the pri- 
mary component in the near-IR using the N90C10 dichroic for 
the observations in the J,H,Ks bands and the JHK dichroic for 
the L, NB4.05 band observations. For sky computation and to 



^ http://www.stsci.edu/software/tinytim/tinytim.html 

^ http://www.stsci.edu/hst/nicmos/performance/photometry/postncs_keywords.] 
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Table 1. Observation log of high-angular resolution imaging of 
the resolved binary Kelu-1 AB 



Date 


Telescope/Instrument 


Filter 


Exp. time 
[sec] 


31/07/2005 


HST/NICl 


F108N 


2560 






CI 1 QM 

r 1 1 jJN 




28/04/2006 


VLT/NACO 


J 


8x80 






H 


8x60 






Ks 


6x40 


26/05/2006 


VLT/NACO 


L' 


8x43.75 






NB4.05 


8x60 


12/08/2006 


VLT/NACO 


Ks 


8x55 


22/02/2007 


VLT/NACO 


Ks 


8x22 


15/05/2007 


VLT/NACO 


Ks 


8x30 


13/03/2008 


VLT/NACO 


Ks 


8x30 


09/05/2008 


VLT/NACO 


Ks 


8x20 


03/06/2008 


VLT/NACO 


Ks 


8x30 



optimize bad pixel/cosmic ray rejection the observations were 
executed in a jitter pattern with 6-8 points (see Table[T]). To per- 
form accurate photometry and astrometry, we chose the nearby 
point source 2MASS J13054388-25403805 from the GSC-II as 
PSF calibrator, which was observed in the same nights with the 
same instrumental settings as the binary. 

On average the observations were executed under good 
atmospheric conditions with a very good seeing of ~ 0.6'' 
FWHM leading to a Strehl ratio typically between 25-35 % 
in Ks. Only the observation on August 12, 2006 had slightly 
worse conditions with a seeing of 0.8'' FWHM and a Strehl ra- 
tio of only 15%, resulting in a larger error estimation for the 
astrometry. 

For the NACO/VLT observations the standard image pro- 
cessing including flat-fielding, bad pixel cleaning, dark and sky 
correction as well as average combination, was accomplished 
using the recommended Eclipse jitter ( |Devillard||1997| ) soft- 
ware package. We obtained very good final images in almost 
all filters with the exception of the images in the NB4.05 fil- 
ter, where we could only barely detect the Kelu-1 AB system 
due to a very low S/N ratio. Thus no photometry was possible 
for this filter. Figure[T] shows the HST as well as the best final 
VLT/NACO images, clearly revealing the increasing separation 
of the binary components with time. 



2.2.2. Spectroscopy 

Our goal to obtain first spatially-resolved spectra of each com- 
ponent of the binary system was achieved by performing spec- 
troscopy with NACO and its H- and ^-filters (spectral range 
1.37 - 1.72yum and 2.02 - 2.53 yum) on April 26, 2006. We used 
the S27_3_SH and S27_3_SK modes, respectively, with a slit 
width of 0.086'' and a resolution of R -1500 and R -1400. A 
total of 20 spectra (in //-band) and 30 spectra (in /T^-band) 
were obtained in an "abbaabb" nodding pattern along the slit 
with an exposure time of 60 s for each frame. The slit was ori- 
ented along the Kelu-1 AB axis, yielding simultaneous spectra 
of both components at each position. As a telluric calibrator 



the B4IV standard HIP 067786 was observed immediately af- 
ter Kelu-1 AB at a slightly lower airmass. 

In a first step the raw, two-dimensional science spectra were 
pairwise subtracted to remove the sky and afterwards divided 
by the median-combined flat-field image. Due to the overlap- 
ping wings of the barely resolved spectra, we used a custom 
made IDL routine to perform an iterative fit to extract the in- 
dividual spectra. This minimum least square fit method turned 
out to be more robust than standard IRAF reduction procedures 
leading to an improved S/N ratio. The wavelength calibration 
was determined based on argon lamp spectra. 

To remove the telluric absorption lines and the instrument's 
relative spectral response function, we used the observed spec- 
tra of the B4 IV telluric standard. After correcting the standard 
star spectrum for its own hydrogen absorption lines, we divided 
it by a blackbody spectrum of Teff = 17,050 K in order to pre- 
serve the continuum shape of Kelu-1 A and B. However, the 
correction with the standard spectra did not sufficiently cor- 
rect for all instrumental artifacts. Therefore we calculated an 
additional response function to correct the remaining slope. 
The added, individual spectra of Kelu-1 A and Kelu-1 B were 
first divided by the unresolved near-IR spectra of Kelu-1 from 
the NIRSPEC Brown Dwarf Spectroscopic Survey (BDSS) 
(McLean et al. 2003| ), which provided a similar resolution. 
Then the single Kelu-1 spectra were divided with this de- 
rived response function, receiving the finally reduced individ- 
ual spectra of Kelu-1 A and B. 

In general the achievable S/N ratio is biased by slit losses 
which are less in H- than in /T-band (wavelength dependent 
since the FWHM of the AO PSF decreases with increasing 
wavelength). The K-hand observation of Kelu-1 AB was in ad- 
dition aff'ected by an alignment error of the binary into the slit, 
which further decreased the S/N ratio. Thus our finally derived 
S/N ratio in the //-band is moderate, though in the /T-band it is 
relatively low, especially for the B -component. 

2.3. Astrometry 

To obtain relative photometry and astrometry of the binary we 
used our custom-made IDL-based simultaneous PSF-fitting all- 
gorithm (B ouy et al.,2003^ , adapted to our HST/NICMOS and 
VLT/NACO dataset. 

In order to recover the binary parameters from the first 
epoch of HST/NICMOS observations in 1998, we tried to fit 
the observed light distribution with theoretical PSFs created 
with Tiny Tim, as well as with observed PSFs. A data base 
of observed NICl PSF stars was compiled by querying the 
HST archive for all observations with NICl in Fl lOM obtained 
within a period of 3 months before or after the Kelu-1 observa- 
tions. To select the best-matching PSFs, only observations with 
NICMOS focus setting "Camera 1", i.e., a PAM (pupil align- 
ment mechanism) focus value close to the value of the Kelu-1 
observations were used. This led to a PSF library of 36 indi- 
vidual NICl exposures in FllOM. Most of the PSFs resulted 
from the observations of very low mass stars and brown dwarfs 
in the Pleiades (GO 7952, [Martin et al.|2000j , which were well 
exposed and represent relatively isolated stars. 



N HST/NIC1 F145M 



.J 



0.2" 
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Fig. 1. Images of Kelu-1 AB obtained with HST/NICMOS and VLT/NACO at different epochs. The orientation and scale are the 
same for all images. Kelu-1 A is the northern component while Kelu-1 B is the component to the south west. 



As the photometry obtained from the spatially resolved NACO 
observations indicates that both components of Kelu-1 are not 
variable by more than a few millimagnitudes, we chose to 
keep the brightness ratio between the two components fixed at 
q = 0.47 for subsequent PSF-fitting, based on the values com- 
puted from the NICl observations in F108N and F113N ob- 
tained in the later epoch. 

For the later obtained HST images of the clearly resolved 
binary, we used a library of 6 diff'erent PSFs for a better error 
estimate: two natural PSFs obtained during previous observa- 
tions of similar objects of the same program and 4 theoretical 
PSFs considering diff'erent focus settings. The relative photom- 
etry and astrometry were measured separately for each of the 
4 images per filter. Finally the results were averaged and the 
uncertainties were calculated from the standard deviation. 
By contrast, for the VLT data only one reference PSF star was 
observed in the same night with the same instrumental settings. 
Therefore all VLT results were based on fitting one single PSF . 
For a better estimate of the accuracy of the results, we repeated 
the PSF fitting using also the previous PSF star observations. 
The derived diff'erences are represented in the error calcula- 
tions. 



2.4. SPITZER Photometry 

The imaging observations with the Spitzer Space Telescope 
(SPITZER) were retrieved from the SPITZER Science 
Archive. As part of the AOR ID 35 (PI G.Fazio) and AOR ID 
3736 (PI B. Goldman) programs, unresolved Kelu-1 AB imag- 
ing data had been obtained on Jan 19, 2004 and Feb 05, 2006, 
respectively, using the Infrared Array Camera (IRAC). IRAC 
provides simultaneous broad band images at 3.6, 4.5, 5.8 and 
8.0 fim (also referred to as channels 1-4 respectively) with a 
pixel scale of 1.2'' and a FoV of 5.2' x 5.2'. Using a dichroic 
beam splitter, two adjacent fields are imaged in pairs (channels 
1 and 3; channels 2 and 4). Then the telescope was nodded to 
image the target in all four channels. 

The SPITZER/IRAC data analysis is based on Spitzer 
Science Center (SSC) pipeline-created mosaic images (with 
pipeline versions S13.2 and S14.0, respectively), the so-called 
post-BCD data (pbdc). The photometry was extracted using 
again the IRAF phot routine. As a first step we converted the 
physical unit of the images from surface brightness into counts 
by dividing with the conversion factor (found in the header of 



each image) and multiplying by the eff'ective integration time of 
268s. For the aperture photometry we chose an aperture size on 
source of 7 pixels and a sky-radius of 7-14 pixels. This guaran- 
teed to receive all PSF flux and still have a sky area without any 
object or bad pixel in it. Since the absolute IRAC calibration 
is based on an aperture with a radius of 10 pixels we applied 
an aperture correction as listed in the IRAC Data Handbook 
3.0. For re-transformation of the net source counts into flux we 
multiplied with the flux conversion factor and the solid angle 
for each pixel. Eventually, the instrumental magnitudes were 
converted into the Vega magnitude scale. The final results are 
based on the combination of the data sets of both programs and 
are listed in Table 13. 3 1 



3. Results and Discussions 

3.1. Orbital parameter determination and dynamical 
mass estimates 

Being clearly resolved as a binary on March 04, 2005 for 
the first time with a separation of 284.3 ±0.8 mas and a po- 
sition angle of 220.9°+ 0.3° ( |Gelino et aL][2Q06 ), the separa- 
tion between the two components increased relatively steady 
to 366.0 ±0.6 mas until June 03, 2008. But since the PA only 
changed by 1.84° during that time, a reliable new orbit predic- 
tion would have been deficient. 

Therefore, we also reexamined the HST/NICMOS obser- 



vations obtained on Aug 14, 1998 (see Sect. 2.1 above) and 
derived a separation of 69.4 ±5.1 mas and a PA = 22.7°± 8.2° 
as best fitting binary parameters. As a separation of 65 mas 
(or less) corresponds to 70% of the diff'raction limit of HST at 
1 . 1 yum, we did not attempt to fit the binary parameters at still 
longer wavelengths. At an observing wavelength of 1.5 yum, 
65 mas are almost exactly half the diff'raction limit of HST. 
Table [2] lists all astrometric results for the Kelu-1 AB system 
obtained so far. 

Compared to the very first orbital predictions by ILiu & 
Leggett| ( [2005 1 ) and Gelino et al. ( 20Q6| ), our nine additional 
data points allow us to improve the determination of the or- 
bital parameters substantially. Based on all astrometric mea- 
surements obtained so far, the best orbital solution was deter- 



mined with our own IDL-based algorithm (see [Kohler et al. 
( 2008 ) for a detailed description of the method). The procedure 
consists of two steps: First, a grid search in period P, eccen- 
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Date 


MJD 


Tele s c op e/In s trunient 


Filter 


Sep. 


FA 


References 










[ nias ] 


rdesl 




14/08/1998 


51039.8 


HST/NICl 


FllOM 


69.4 + 5.1 


22.7 + 8.2 


this work 


04/03/2005 


53434.0 


Keck/NIRC2 




284.3 + 0.8 


220.90 + 0.30 


"Gelino et al."( 20061 


01/05/2005 


53491.0 


Keck/NIRC2 




291.0 + 2.0 


221.20 + 0.60 


Liu & Leggett (26^ 


31/07/2005 


53582.0 


HST/NICl 


F108N 


299.8 + 0.2 


221.33 + 0.04 


this work 


28/04/2006 


53853.1 


VLT/NACO 


Ks 


320.3 + 1.7 


221.84 + 0.1 


this work 


26/05/2006 


53881.1 


VLT/NACO 




322.7 + 1.9 


221.86 + 0.1 


this work 


12/08/2006 


53959.9 


VLT/NACO 


Ks 


330.6 + 5.0 


222.24 + 0.1 


this work 


22/02/2007 


54153.4 


VLT/NACO 


Ks 


341.5 + 1.5 


222.39 + 0.1 


this work 


15/05/2007 


54235.9 


VLT/NACO 


Ks 


345.3 + 1.6 


222.48 + 0.1 


this work 


13/03/2008 


54538.8 


VLT/NACO 


Ks 


362.8 + 0.6 


222.63 + 0.1 


this work 


09/05/2008 


54595.7 


VLT/NACO 


Ks 


364.5 + 0.5 


222.67 + 0.1 


this work 


03/06/2008 


54621.5 


VLT/NACO 


Ks 


366.0 + 0.6 


222.74 + 0.1 


this work 



of 177^^^^ Mjup. This is the first purely dynamical mass esti- 
mate for the Kelu-1 AB system with no evolutionary models 
involved in the calculations. 

This mass prediction reveals a slightly higher total sys- 
tem mass than one would expect for a BD-BD binary system, 
where an upper mass limit of 72 - 75 Mjup for brown dwarfs 
(the substellar mass-boundary) would result in a maximum to- 
tal system mass of -150 Mjup. Though, within the uncertain- 
ties of our estimate a pure BD-BD binary is still possible. 
The known presence of Li absorption in the unresolved opti- 
cal spectrum (Ruiz et al.||1997[ Kirkpatrick et al.|[T999j indi- 
cates that at least one component must have a mass ^ 65 Mj^p 
which is the lithium-burning limit ( Rebolo et al.||1992{ Basri 
1998| ). If the system consists of two objects, this would imply 
that the primary has a mass of ^ 1 10 Mjup, which in turn would 
correspond to a main-sequence spectral type of M8V. This is 
in contradiction to its early L spectral type as it has been de- 
termined both from integrated spectra ( [Kirkpatrick et al.|1999[ 
Dahn et aT"2002) as well as in §|3.2| below from the spatially 



tricity e and time of periastron Tq is carried out, while Singular 
Value Decomposition is used to solve the linear equation sys- 
tem for the four Thiele-Innes constants. From the Thiele-Innes 
constants, the remaining orbital elements (semi-major axis a, 
argument of periastron oj, position angle of the line of nodes fl, 
and inclination /) are computed. The result of the grid search is 
as function of P and e, which helps to assess the quality of 
the fit, and the confidence limits of the fit parameters. 

In the second step, all 7 orbital elements are fit simultane- 
ously by a Levenberg-Marquardt algorithm. This further im- 
proves the fit by interpolating between grid points and ensures 
that all 7 elements are treated equally in the final fit. To avoid 
that the algorithm converges on a local instead of the global 
minimum, we decided to use all orbits resulting from the grid- 
search as starting points and carried out 148 x 100 runs of 
the Levenberg-Marquardt algorithm (the grid spans 148 peri- 
ods and 100 eccentricities). 

To derive estimates of the confidence intervals for the pa- 
rameters, we studied the function around its minimum. To 
find the confidence interval for one parameter (for example Tq), 
we perturbed this parameter away from the minimum in;^^ and 
optimized all the other parameters. Any perturbation of a pa- 
rameter will of course lead to a larger x^- The range in Tq 

within which x^(To) - xl^n < ^ defines the 68 % confidence Table 3. Parameters for the best orbital solution of Kelu-1 AB 
interval for Tq. This interval is usually not symmetric around 
the To of the best fit, therefore we list in Table [3] separate limits 
for positive and negative perturbations. 

Finally, the total system mass is computed from the semi- 
major axis and orbital period through Kepler's third law. To 
convert the semi-major axis from arcsec into AU as well as to 
compute the mass in Mjup we used a distance of 18.7 ± 0.7 pc 
as determined by |Dahn et aL] ( |2QQ2| . 

We applied this IDL code to the relative astrometric data 
as given in Table [2] (hereafter data set N°l) and found a global 
minimum reduced of 1.3. The best-fit results are summa- 
rized in Table |3j including the uncertainties of each fitted pa- 
rameter, and the orbital solution is illustrated in Fig. [2] The 
orbit is quite eccentric (e = 0.82 ± 0.10 ) with an inclination of 
84.9 



resolved spectroscopic observations of Kelu-1 A and B. 



Parameter 


Value 


System mass Ms [Mjup] 


177 +1^3 

^ ^ -55 


Period P [years] 




Semi-major axis a [mas] 


339 :>f 


[AU] 




Eccentricity e 


0.82 


Argument of periastron oj [°] 


C7 o +15.0 
^ -20.0 


PA. of ascending node Q. [°] 


39 4 +2-^ 


Inclination / [°] 


84.9 


Date of periastron Tq 


2452079 ^200 



hl.O 



2Q , i.e. seen almost edge-on. The period of 38^^ years 
and the semi-major axis of 6.4^ ^ ^ AU yield a total system mass 



reduced 



(Jun 18, 2001) 
1.3 
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offset in RA [mas] 

Fig. 2. Orbit of the Kelu-1 AB system. The star indicates the position of Kelu-1 A while the measurements for the relative motion 
of Kelu-1 B around Kelu-1 A are indicated by crosses. Measurement errors are comparable to or smaller than the plotted symbols. 
The resulting best orbit fit is represented by the solid, blue curve and the green, dotted lines show the orbits at the 90% confidence 
limit (with periods of 32 and 46 years). The red, dash-dotted line marks the line of nodes and the dashed line indicates the position 
of the periastron. In addition the positions on January for the years 1998-2008, as well as the predicted position for 2012 are 
marked. 



3.1 .1 . Comparison to previous orbit predictions 

Since the HST/NICMOS observation from 1998 is obviously a 
key component to fit the orbit, we compared our binary param- 
eters for this observation with those obtained by [Gelino et al.| 
( 2006| ). They derived a separation of 45 ± 18 mas and a PA = 
38.0° or 21 8.0° ± 1 1 .9° (due to their 1 80° ambiguity in position 
angle due to the unresolved objects). Within the uncertainties, 
these values are in good agreement with our measurements. In 
an additional check we analyzed what orbit solutions we would 
get with their slightly diff'erent result. We applied our IDL code 
on two additional data sets which both contained the same val- 
ues as data set N^l except for the very first measurement. In 
data set N°2 we replaced the values of the August 1998 obser- 
vation by a separation of 45 ± 18 mas and a PA of 38.0° while 
data set N°3 used a PA of 218.0° . The best-fitting orbital solu- 
tion for data set N°3 results in a physically unacceptable total 
system mass of 108 solar masses. If we consider only orbital 
models with a total system mass ^ 230 Mjup (assuming a possi- 
ble triple system), we find a minimum reduced of 1.5, which 
is higher than that for data set N° 1 . Therefore we exclude the 
PA of 218° as a possible first data point. The result from the 
best-fit solution for set N°2 yields, without any constraints, a 
system mass of 97^23^ Mjup, period of 28^^^ years and semi- 
major axis of 4.3^Q 2 AU. These values are all slightly smaller, 
though within the uncertainties still consistent with the result 
for data set N° 1 . However, since its minimum reduced of 
1 .4 is still higher than the one for set N° 1 we stay with our own 



astrometric result for the first HST/NICMOS data point. In ad- 
dition, this ensures a consistent data reduction method through- 
out the whole monitoring program. 

Comparing the final results in Table [2| with the previo usly 
published orbit parameter predictions by jLiu & Leggett| and 
[Gelino et al . we first exclude their assumed circular orbit, since 
with an eccentricity < 0.01 this would result in unphysical pa- 
rameters like a dynamical total system mass of 1.06 M© or 
a reduced > 3. Nevertheless, our derived inclination of 
84.9^^-^° agrees well with the limiting values of 81 < / < 90° 
reported by [Gelino et al.[ With a total system mass of 177 Mjup, 
our prediction is higher than the mass estimates of ~ 120 Mjup 
and -115 Mjup by Liu & Leggett and Gelino et al. respectively. 
However, their calculations were purely based on photometric 
results compared with theoretical models and a still uncertain 
age estimation for the system. 

Considering a possible third component in the system (as 
derived from the results in § [3.2.1| below), would allow all ob- 
jects to have component masses < 95 Mjup, and hence spec- 
tral types of LO or later. This would resolve the above men- 
tioned discrepancy between a system mass of -180 Mj^p, spec- 
tral types of LO and later, and the fact that at least one of the 
components must have a mass less than 65 Mj^p as indicated 
by the presence of Lithium in the combined spectra. 

All together we showed a possible solution in reasonable 
agreement with the observations, but most likely not the final 
one, since the fraction of the orbit observed so far is too small 
for a reliable and consistent orbit determination. For a robust 
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Fig. 3. Left: H- and ^-band spectrum of Kelu-1 A (smoothed to template resolution, solid line) compared to the best-fit binary 
composites of 2M2107-0307 (LO) combined with 2M0727-1710 (T7) (red, solid) and 2M1439-1929 (LI) with 2M0415-0935 
(T8) (blue, dotted). In addition the single object spectrum of 2M2 107-0307 (LO) (green, dashed) is plotted for comparison. Right: 
close-up of the //-band and the peculiar feature starting at 1.6 jim. While the single LO (blue, dotted) and LI (green, dashed) show 
a plateau shape from 1.55 - 1.64 jim due to the FeH absorption band heads, the LO -h T7 binary composite (red, solid) reproduces 
the deep depression to a much higher accuracy. 



orbital solution, more observations are required, at least up to 
the turning point in the (projected) orbit expected around the 
year 2012 (see Fig.|2]). Therefore further investigations, includ- 
ing the measurement of the radial velocities of the system com- 
ponents, will be done with our ongoing long-term monitoring 
program. 

3.2. Spectral types 

The extracted and reduced HK- spectra of Kelu-1 A and Kelu- 
1 B are shown in Figs. |3] and |4] respectively. To determine the 
spectral type of each component we compared them with low- 
resolution spectra of 5 late M dwarfs (M8-M9.5) and 39 field L 
dwarfs (L0-L9) from the SpeX prism data archiv^ All spectra 
were normalized to the //-band. In addition the Kelu-1 spectra 
were smoothed to the template resolution and scaled by a factor 
of 0.95 to eliminate the normalization bias. 

3.2.1. Kelu-1 A 

The general //-band slope of Kelu-1 A best fits between the 
optically classified LO 2MASSI J21073 16-030733 (hereafter 
2M2107-0307, ,Burgasser et al.,2004j and the LI optical stan- 
dard 2MASSW J14392 83-F192915 (hereafter 2M 1439-^1929, 
[Burgasser et al.|[2004| ). In the /T-band the overall shape is 
clearly very similar to that of 2M1439-i-1929 with the same 
strong CO absorption band heads at 2.3-2.4yum and, at the 
short- wavelength end of the ^-band, a steeper slope than in the 
LO spectrum, due to stronger H2O absorption. This suggests a 
spectral type of L0.5 ± 0.5 which is slightly earlier than the 
derived spectral type of LI. 5 - L3 from Liu & Leggett ( 2005| ) 



http://www.browndwarfs.org/spexprism 



whose estimates were only based on JHK colors and absolute 
magnitudes. However, the //^-band spectrum of Kelu-1 A shows 
a pronounced feature, a deep, broad depression from 1.6 - 
1.65 yum. This feature is slightly off'set from the 1.55 - 1.64yum 
FeH absorption band heads which produce a flat plateau in this 
wavelength range and are usually present in early to mid type 
L dwarfs ( Gushing et al.|[20Q3] ). Instead, the 1.6yum depres- 
sion corresponds much more to the CH4 absorption band which 
starts to appear at the L/T transition and strengthens through the 
T sequence and is therefore very unusual for an early L dwarf. 
In addition there is no indication of the 2.2yum CH4 absorption 
line, another defining feature for the T dwarfs. 

Recently the same feature has been observed in 3 
formerly classified single, early to mid-type L dwarfs 
(2MASS J05 185995-2 828372, |Cmz et al.| [20041 SDSS 
J08053 1 . 84-^48 1233.0, |B urgasser'2007b'; 2MASS J03202839- 
0446358, Burgasser et al. 2008 ) . By combining an early to mid 
L dwarf spectrum with an early to mid T dwarf spectrum, these 
objects were reinterpreted to be unresolved spectroscopic bina- 
ries, with the depression being very well reproduced by a com- 
bination of the FeH absorption in the L dwarf primary and the 
CH4 absorption of the T dwarf secondary. Two of the targets 
have been resolved by observations since then: 2MASS 0518- 
2828 with high-resolution imaging ( [Burgasser et al.|2006b| ) and 
2MASS 0320-0446 with high-resolution spectroscopy to be a 
single-lined spectroscopic binary (SBl) ( [Blake et al.|2008[ ). 

To probe if we can also reproduce the depression in the 
Kelu-1 A spectrum including the overall shape of the H- and 
/T-band, we created a large set of binary spectral templates by 
combining late M-dwarf / early L dwarf spectra (M8-L4) with 
early to late T dwarf spectra (T2-T8) from the SpeX prism 
data archive. We selected the SpeX data archive since it is the 
only archive which provides a continuous and large set of spec- 
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Fig. 4. Left: H- and ^-band spectrum of Kelu-1 B (smoothed to template resolution, solid line) compared to the best-fit spectra 
of the red L dwarfs 2M0141-4633 (LOpec, red solid line) and 2M0 103-^1935 (L6, blue dotted line), as well as the field L dwarf 
2M1506-1321 (L3, green dashed line). Right: close-up of the //-band with the same objects as on the left. The distinct FeH 
absorption dip around 1.58yum is present in the L3 spectrum, but not in the LOpec or L6 spectra, while the overall slope better 
fits the red objects. 



tral types from late M-dwarfs to the latest T dwarfs in the H- 
and /T-band and, therefore, provides the most homogenous set 
for assembling all the diff'erent binary templates. After nor- 
malization we calculated the RMS between the science spec- 
trum and the binary templates to determine the best matches. 
The results with the least deviation are displayed in Fig. [3] 
(left). The combination of 2M2107-0307 (LO) and 2MASS 
J07271824+ 17 10012 (T7 standard) (2M0727-hl710, .Burgasser^ 
et al.|2006a| best fits the slope and depression in the H-hmid. 
Whereas the smallest RMS in the K-band is achieved with 
the combined spectrum of 2M1439-1929 (LI) and 2MASS 
J0415 1954-0935066 (T8 standard) (2M04 15-0935, [Burgasser 



et al. 2004[ ) with a better match of the steeper slope at the short- 
wavelength end of the /T-band spectrum due to stronger H2O 
absorption. On the right side of Fig. |3]we show a close-up of 
the //-band and the peculiar feature starting at 1.59yum. One 
can clearly see that neither the single best-fit LO nor the LI 
match the deep depression while the L0-hT7 binary template re- 
produces the whole spectra to a much higher accuracy. Beside 
those best-fit comparisons there were several other matches 
with similar good results among the binary templates with sec- 
ondaries of type T6.5 and later. Including the small diff'erence 
between the best match with an LO primary in //-band and an 
LI primary in /T-band we derive a component spectral type of 
L0.5 ± 0.5 for Kelu-1 Aa and T7.5 ± 1 for Kelu-1 Ab. 



3.2.2. Kelu-1 B 

A precise spectroscopic classification of Kelu-1 B is even 
more challenging. The object shows an //-band continuum 
notably diff'erent from the usual plateau- shape of "normal" 
early to mid- type L dwarfs (see Fig. |4]). Its sharply peaked, 
triangular-shaped continuum resembles much more that of 
the young, red LOpec 2MASS J01415823-4633574 (2M0141- 



4633, [Kirkpatrick et al.| [20061 ) and the red L6 2MASSI 
J0103320-F193536 (2M0103-F1935, |Cruz et al.1[2Q04l ). While 
the distinct dip around 1.58yum in the Kelu-1 B spectrum 
(caused by two prominent FeH absorption bands) is not visi- 
ble in the red spectrum of the LOpec and L6 dwarf, it is in turn 
present only slightly weaker in "normal" L dwarf spectra. 

By contrast the /T-band spectrum of Kelu-1 B does not fit 
any of the two red objects. While its overall flux distribu- 
tion fits relatively well that of the L3 optical standard 2MASS 
J15065441-hl321060 (2M1506-hl321, |Burgasser| |20Q7a|, it 



again reveals a triangular- shaped continuum with a slightly red- 
ward shifted /T-band flux. In addition the CO absorption bands 
are somewhat weaker as in the L3 dwarf. 

For comparison in the //-band, only the above mentioned 
two unusual red L dwarfs were available in the archive. 
Therefore, we used the surface gravity independent H2O 
indices in the short- wavelength end of the //-band (1.49- 
1.56 yum) from Aller s et aL] ( |2007| ) for a better estimation of the 
spectral type. In a first attempt we applied this H2O index to 
the LOpec 2M0104-4633 and obtained a spectral type of M9.7 
± 0.5, which is fully consistent with the LO classification of 
[Kirkpatrick et al.|(2006| ). Using the same indices for Kelu-1 B 
we receive a spectral type of L2.6 ± 1.5 where the uncertainty 
includes the scatter of the index relation and the error of our 
index measurements. Using the results from our fits and calcu- 
lations we thus adopt an overall spectral type of L3pec ±1.5 
for Kelu-1 B. 

Similar peaked //-band spectra and unusual red near-IR 
colors have been so far mainly detected as characteristics for 
young ultracool dwarfs ( [Luhman et al.|2004 Kirkpatric k et al 



2006 i AUers et al.|20071 [McElwain et al.|20b7 t Metchev et al 



[2008|. The eff'ect can be explained by two possible scenarios 
which are both caused by lower gravity and lower atmospheric 
pressure: (i) a diminishing collision-induced absorption (CIA) 
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of H2 in the H- and TT-band ( [Kirkpatrick et al.||20Q6l ), or (ii) 
increased water vapor absorption on each side of the H-band 
peak (ILuhm an et al.|2004t[Allers et al.|20Q7| ). Another possible 
explanation for the peculiar shape is a higher metallicity which 
implies a higher opacity caused by greater dust production. 
This again results in an overall redder spectrum with weaker 
alkali (Na, Ki ) and molecular (FeH and H2O) absorption lines 
and a red-shifted ^-band peak. 

While Kelu-1 B also shows slightly redder J-Ks colors (see 
§ |3.3| below), the triangular-shaped H- and /sT-band spectra can 
not be explained by low gravity and thus implying a young age, 
since Kelu-1 A, the other member of the coeval binary system, 
shows (even with a possible companion) the "normal" behavior 
of an early L dwarf and no spectral evidence of a young ultra- 
cool object. Recently |Looper et al.| (|2Q08|) identified two new 



peculiar L dwarfs which showed the same triangular- shaped H- 
and ^-band spectra as Kelu-1 B. For one of those objects, the 
L6.5 pec 2MASS J21481628-h4003593 they suggest that the ef- 
fect may arise from a metal rich atmosphere rather than from 
low gravity, since its large tangential velocity does not indicate 
an unusual young age. 

Since low gravity as well as higher metallicity result in 
weaker alkali and FeH absorption, especially in the red-optical 
and J-band spectra, further investigations on resolved Kelu-1 B 
spectra in these wavelength ranges will be necessary for a bet- 
ter spectral characterization. 



3.3. Magnitudes and Colors 

We calculated the individual magnitudes of the components in 
the J, H, Ks and U filters from the measured binary flux ra- 
tios (see § [23]) and the integra ted 2MASS ( jCutri et al.|(2QQ3] ) 
and UKIRT ( Leggett et al .^2002 ) magnitudes of the unresolved 
Kelu-1. Even if the 2MASS and VLT/NACO near-IR filters are 
not exactly identical we did not apply any correction factor, 
given that the spectral energy distributions (SEDs) of the Kelu- 
1 AB components are so similar that the flux ratio should not 
be noticeably afl'ected by the difl'erence in these two photomet- 
ric systems. For the UKIRT U filter one has to consider that it 
cuts on slightly bluer wavelength than the VLT/NACO filter. 
However, [Leggett et al.| pOQ2 ) showed that a similar diff'erence 
between the UKIRT and MKO filter in this particular wave- 
length regime is only significant for late L dwarfs and T dwarfs 
when the CH4 absorption at 3.3 yum starts to play a role. For 
early L dwarfs they claim the eff'ect to be less than the uncer- 
tainties in the published magnitudes. Thus we used no correc- 
tion factors between the UKIRT U and NACO U filter system. 

Table |4] summarizes the photometric results. For the de- 
termination of the absolute magnitudes we used a distance of 



18.7 + 0.7 pc from Dahn et al. ( 2002 ). Our error estimations 
include the uncertainties in the flux ratios, the corresponding 
uncertainties of the "unresolved" magnitudes, as well as the 
error of the distance. For comparisons with J, H and K magni- 
tudes and spectral type relations derived in the MKO photomet- 
ric system, we used the transformation equation for 2MASS 




Fig. 5. Kelu-1 AB image achieved with VLT/NACO on 26. 
May 2006. This is the first resolved L'-band image of a brown 
dwarf binary system. 



The photometric monitoring in the Ks-band over the past 
2 years revealed no evidence of photometric variability within 
the errors for Kelu-1 A or Kelu-1 B. Compared to the values 



of Liu & Leggett ( 2005| , our individual magnitudes in H and 
Ks-band are in good agreement within the errors, whereas our 
J-band magnitude for Kelu-1 B is slightly fainter. This is likely 
due to in general lower AO performance in J-band, resulting in 
a considerably worse Strehl ratio and leading to broader PSF 
wings, which were hence not separated enough for an accurate 
flux ratio calculation. We took this into account by assuming a 
larger error. While the resulting J-Ks and J-H colors for both 
components are at the red end, both H-Ks colors fit very well 
into the color range for spectral types of LO-Ll and L3-L4 (see 
[Cushing et al.||2008| ) respectively. The red colors of Kelu-1 B 
correspond to the redder H- and ^-band spectra reported and 
discussed in § |3.2.2| 

We note that our derived L'-band photometry with 
L\ = 11.42 ±0.15 mag and L'b = 11.66 ± 0.15 mag are the first 
resolved L-band magnitudes for any L dwarf binary so far (see 
Fig. |5]). While the unresolved U photometry was as overlu- 
minous as in the other near-IR bandpasses for a spectral type 
of L3 ( [Golimowski et al.||2004b| , our resolved (Kmkq-LOa 
= 0.88 + 0.15 and (Kmko-LOb = 1.03 + 0.15 colors, as wefl 
as the absolute U magnitudes are consistent with a spectral 
type of L0.5 and L3, respectively (as in § |3.2| ), derived from 
the color-spectral type and magnitude- spectral type relations in 
[Golimowski et al.| ( |2QQ4b| ). 

To get an idea of the photometric values of the potential 
third component in the Kelu-1 AB system, we computed the J, 
H and Ks magnitudes from the Kelu-1 A magnitudes. We con- 
sidered the allowed flux ratios of Aa and Ab in a filter, as well 



as the latest 2MASS color- spectral type relations in Cushing 
et al.| p008 ) under the assumption of spectral types ~ L0.5 and 



^ T7 as derived in 



3.2[ The primary Kelu-1 Aa would be only 

13.89, 



into MKO of Stephens & Leggett ( 2004 ). 



slightly fainter than the composite source with Ja 
HAa= 12.95 and (Ks)Aa = 12.33 while Kelu-1 Ab is notably 
fainter: JAb = 16.99, HAb = 16.95 and (Ks)Ab = 17.07mag. We 
do not assign any uncertainties to our magnitudes, since there 
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Table 4. Apparent and absolute magnitudes of Kelu-1 AB 



Telescope/Instrument 


Filter 


Kelu-1 AB 

111 Liiiagj 


A mag 


Kelu-1 A 
m [mag] M [mag] 


Kelu-1 B 
m [mag] M [mag] 


HST/NICl 


F108N 




0.83 + 0.02 


14.56 + 0.02 


13.20 + 0.08 


15.31+0.02 


13.95 + 0.08 




F113N 




82 + 02 


14.28 + 0.02 


12.92 + 0.08 


15.01+0.03 


13.65 + 0.09 


2MASS ^ 


J 


13.414 + 0.026 












VLT/NACO 


J 




0.82 + 0.09 


13.83 + 0.08 


12.47 +0.11 


14.65 + 0.12 


13.29 + 0.14 


2MASS ^ 


H 


12.392 + 0.025 












VLT/NACO 


H 




0.51 +0.02 


12.92+ 0.03 


11.56 + 0.09 


13.44 + 0.04 


12.08 + 0.09 


2MASS ^ 


Ks 


11.747 + 0.023 












VLT/NACO 


Ks 




0.39 + 0.01 


12.32 + 0.02 


10.96 + 0.08 


12.72 + 0.03 


11.36 + 0.09 


UKIRT/IRCAM ^ 


U 


10.78 + 0.15 












VLT/NACO 






0.24 + 0.01 


11.42 + 0.15 


10.06 + 0.17 


11.66 + 0.15 


10.30 + 0.17 


Spitzer/IRAC 


3.6 fim 
4.5 jim 
5.8 yL/m 
8.0 yum 


10.90 + 0.03 
10.85 + 0.04 
10.68 + 0.07 
10.55 + 0.05 













^ from 
^ from 



Cutri et al.|([2003]) 



Leggettet al. (2002i 



was no direct photometric measurement and we just want to 
give a first prediction of possible photometric values. In addi- 
tion the results are in a very good agreement to the absolute 
magnitude- spectral type relations of |Liu et J^ ( |2QQ6| ). 



3.4. Testing evolutionary tracks 
3.4.1. Ages 

Age estimates of field brown dwarfs without established clus- 
ter membership have to rely on theoretical isochrones, using 
the mass-luminosity-age relation for brown dwarfs. Figure [6] 
displays the color-magnitude diagram of Kelu-1 A and B com- 
pared to the widely used DUSTY models of Chabrier et al.| 
( |2Q00| ) for early to mid-type L dwarfs. For Kelu-1 A these mod- 
els predict an age clearly between 0.1 and 0.5 Gyr taking into 
account the estimated photometric errors. For Kelu-1 B, how- 
ever, they predict a slightly older age of 0.3 - 1 Gyr. Since bi- 
nary components are believed to be coeval, we assess a min- 
imum system age of 0.3 and a new upper age limit of 0.5 
Gyr. In addition, the triangles in Fig. |6] show the position of 
of the Kelu-1 Aa and Kelu-1 Ab components if they would be 
resolved. Therefore we included the COND models of lBaraffel 
|et al.| ( |2003| ) which better consider the dust condensation of T 
dwarfs. The resulting age for both individual components is 
consistent with the above estimation. The here derived age es- 
timate (0.3-0.5 Gyr) further decreases the upper age limit for 
Kelu-1 AB compared to previous estimates from [Basri et al.| 
([1998 ) (Q.3-lGyr ) and|Liu & Leggett| (|2005]) (03.-0.8 Gyr). The 
determination by |Basri et aL ( 1998| ) was solely based on the 
strength of the unresolved Lii absorption feature, while |Liu"&| 
|Leggett| ( [2005 1 ) relied on the mass evolutionary tracks of the 
same model. They took the presence of lithium and the sub- 
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Fig. 6. Color - Magnitude diagram: The red asterisks show the 
location of the resolved components Kelu-1 A and Kelu-1 B. 
The position of the components Kelu-1 Aa and Ab, if they were 
resolved, are indicated by the dotted lines and the triangles. 
While the DUSTY isochrones are displayed as solid lines in 
the upper right part, the COND isochrones are displayed by the 
dashed lines. 



stellar nature of both components into account, which restricts 
the maximum mass for each component. 

3.4.2. Masses 

In order to compare our dynamical mass determination of 
Kelu-1 AB with mass predictions from evolutionary models, 
we again utilize the DUSTY and COND models as a function 
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Fig. 7. The absolute Ks magnitudes of the Kelu-1 system com- 
ponents compared to the DUSTY (solid lines) and COND 
(dashed lines) models at 0.05, 0.1, 0.5, 1 and 5 Gyr to derive 
maximum mass estimations for each component. 

of absolute Ks magnitude and age (see Fig. [7]). Using the up 
per limit for the system age of 0.5 Gyr as obtained in § 3.4.1 
we derive maximum individual masses of 61 ± 2 Mjup for Kelu- 
1 A and 50 + 2 Mjup for Kelu-1 B resulting in a maximum total 
system mass of 111 ±4Mjup. Even if we include the possible 
third component Kelu-1 Ab, the additional mass of 18.5 Mjup 
would only increase the total mass to ~ ISOMj^p. This mass is 
still about 47 Mj^p lower than the total system mass we derived 
in § |3.1| by dynamical mass determination. This discrepancy 
between luminosity, mass and age in the evolutionary tracks 
might result in an under-prediction of the masses. A similar 
finding has already been reported for other objects. Higher sys- 
tem masses were derived by dynamical mass determination, by 
[Close et al.| (i2005) for AB Dor C, Ireland et al.lf2008) for GJ 
802 b or just recently [ Cardoso eTal . ( 2008 ) and Kasper et al. 
(ApJL, 2008 submitted) for e Ind Ba, Bb. They all claim as pos- 
sible reasons the underestimation of the abundance of conden- 
sates in the atmospheres as well as the eff'ect of magnetic fields 
hindering heat flow or even the unusual assumption of non- 
coevality in binary systems. 

4. Summary and outlook 

We presented new astrometric and photometric observations 
of Kelu-1 AB. For the first time we achieved resolved spec- 
troscopy of the individual components of this brown dwarf bi- 
nary system. The near-IR spectra display interesting peculiar- 
ities for each component. The spectrum of Kelu-1 A shows a 
distinct cup- shaped dip from 1 .6 -1 .65 jim which did not fit any 
early, single L dwarf spectrum, but could be best reproduced by 
a combined spectrum of an early L- and late T dwarf. Our spec- 
tral template matching algorithm leads to the result that Kelu- 
1 A itself is a close spectroscopic binary with spectral types of 
L0.5 ± 0.5 and T7.5 ± 1. In contrast, Kelu-1 B is characterized 
by a triangular- shaped continuum, similar to very young red L 
dwarfs and was classified as L3 pec ± 1.5 . 



Eight epochs of high-resolution astrometry from 2005- 
2008 enabled us to compute the orbital parameters and de- 
rive the first dynamical mass measurement for the system. We 
found a quite eccentric orbit (e = 0.82 + 0.10 ), seen almost 
edge-on. The period of 38^^ years and a semi-major axis of 
6.4^^ 3 AU yields a total system mass of 177^^^^ Mjup which 
is slightly higher than one would expect in a BD-BD binary 
system. If the system consists of only two objects, this would 
imply that the primary is more likely a very low mass star with 
a mass ^ 85 Mj^p. This is in contradiction to its defined early 
L-spectral type from spectroscopy which, along with its age 
estimation, indicates an affiliation to brown dwarfs. Including 
a possible third component in the system would solve this dis- 
crepancy with all component masses being below the substellar 
mass limit of 75 Mjup. The non-detection of this component in 
any previous high-resolution imaging observation (Keck, HST, 
NACO), indicates a maximum separation of ~ 0.22'' 4. 1 AU) 
related to the assumed brightness diff'erence of Kelu-1 Aa and 
Ab. 

The derived magnitudes of L' = 1 1 .42 ± 0. 15 mag for Kelu- 
1 A and L' = 1 1 .66 ± 0. 15 mag for Kelu-1 B provide the first re- 
solved L-band photometry for any L dwarf binary. These values 
and the related near-IR colors support the spectral type deter- 
minations from spectroscopy. Comparing the derived near-IR 
photometry of the individual components with the widely used 
theoretical DUSTY models, led to a slightly younger age esti- 
mation of 0.3 - 0.5 Gyr compared to previous estimations of 
its age. Using this new upper age limit, the evolutionary tracks 
yield only a maximum total system mass of -130 Mj^p, even 
if a third companion is considered. This arising discrepancy to 
the dynamically determined mass shows the importance of fur- 
ther mass determinations independent of theoretical models. 

With all the above conclusions, the hypothesis that Kelu- 
1 A itself is a binary system seems conclusive. If the third com- 
ponent in the Kelu-1 AB system can be confirmed in further 
observations, this would be the first pure brown dwarf triple 
system detected so far. 

Future work includes resolved optical spectra of Kelu- 
1 AB. This will ascertain in which component Li absorption 
is present and thus help to better establish the ages of the com- 
ponents. In addition, it will extend the individual spectral type 
determinations into the red optical. This might give further 
conclusions on mechanisms like low-gravity, high-metallicity 
and/or cloud condensation leading to the redder, triangular 
shaped near-IR spectra of Kelu-1 B. The continued monitor- 
ing of Kelu-1 AB will help to refine the orbital parameters and 
its dynamical mass. Further, upcoming spatially resolved high- 
resolution radial velocity observations will provide information 
on the spectroscopic orbit for the determination of individual 
masses and possibly resolve the third component. These re- 
sults will help to improve and adjust the theory of brown dwarf 
evolution and improve the interpretation of the observations of 
these objects. 
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